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We compare in detail the results of simulations of electromagnetic showers in ice in the GeV- 
TeV energy range, using both the GEANT package and the ZHS Monte Carlo, a code specifically 
designed to calculate coherent Cerenkov radio pulses from electromagnetic showers in dense media. 
The longitudinal and lateral profiles as well as the tracklengths, and excess tracklengths are shown to 
agree at the 10% level. We briefly comment on the negligible influence of the Landau-Pomerancuk- 
Migdal effect on the total shower tracklength. Our results are relevant for experiments exploiting 
the radio Cerenkov technique in dense media as well as for other detectors that rely on the Cerenkov 
effect in water or ice. 

PACS numbers: 96.40.Pq, 95.85.Bh, 95.85.Ry, 29.40.-n 



I. INTRODUCTION 

Ultra high energy neutrino (UHEj/) detection is one 
of the experimental fields in astroparticle physics that 
has received most attention in the last two decades [0. 
Several experiments are already taking data and several 
more are under way or in the proposal stage. Due to the 
low neutrino interaction probability and the low expected 
neutrino fluxes, immense volumes of detector material 
(~ 1 km 3 at least) are required B. Currently explored 
detection techniques exploit the observation of Cerenkov 
radiation in the optical frequency range from neutrino- 
induced showers and neutrino-induced charged leptons 
in dense, transparent media g, ||, [| or the search for 
horizontal air showers Q. 

The observation of coherent Cerenkov pulses in the 
MHz- GHz frequency range from neutrino induced show- 
ers in transparent, dense media, provides an alternative 
method of detecting UHEz/s ||. The technique is most 
promising at neutrino energies >PeV (10 15 eV) at which 
effective volumes in excess of 1 km 3 can be achieved in a 
cost-effective manner [Q, [| . Electromagnetic showers are 
known to develop an excess of negative charge of about 
20% mainly due to Compton scattering. When the wave- 
length of the radiation is larger than the typical dimen- 
sions of the shower the emission from the excess charge 
is coherent, and the power in radio waves scales as the 
square of shower energy as predicted by Askary'an in the 
1960's §. 

As interest in high energy neutrino detection grew in 
the mid 1980's, proposals were made to search for ra- 
dio pulses produced by the showers that develop when 
UHE neutrinos interact. Arrays of antennas could be 
used to detect the pulses produced in deep ice and 
radiotelescopes for those produced under the Moon sur- 
face 0| . Clearly a reliable calculation of radio pulses was 
needed to explore the possibilities of these techniques. A 
fast Monte Carlo code to simulate up to PeV electro- 



magnetic showers was specifically designed in the early 
1990's to calculate the interference pattern in the MHz- 
GHz region from an electromagnetic shower in ice (the 
ZHS code from now on) jlj, [l3) . The results revealed the 
wealth of information that is kept in the radiation pat- 
tern |47j . Perspectives for the technique became most 
encouraging for energies above the 10 PeV scale. The ra- 
dio technique was further studied in the 1990's from three 
perspectives: prospects for radio detection arrangements 
were discussed || , experimental measurements were per- 
formed pH ], and the pulse simulations were extended to 
the most promising EeV range where full simulations 
are out of question with conventional computing facili- 
ties H ||. 

The existence of the charge excess was recently con- 
firmed in a revealing accelerator experiment at SLAC 
]l9[ , giving a good thrust to the technique and generat- 
ing a large number of proposals, some of them already 
in operation pc| . Two experiments are presently active: 
the RICE experiment, an array of antennas buried in 
the transparent polar ice cap |l4| and the GLUE exper- 
iment which uses the visible side of the Moon as target 
for UHEi/ and cosmic ray interactions 21 . New propos- 
als include: ANITA |^(|, a balloon flown antenna looking 
down to the polar ice cap; SalSA p2[ , the same concept 
as RICE but exploiting the excellent optical properties 
of some salt domes; a proposal to analyze data from the 
FORTE satellite || @; and the LOFAR project, which 
will use a radio-astronomy array of low frequency anten- 
nas to search for radio emission from extensive air show- 
ers @. 

Radio patterns have been studied for long by an in- 
dependent group, both from the theoretical p5| and the 
phenomenological sides ||, 26 1. Recently the predictions 



of ZHS have been challenged by shower simulations per- 
formed with the all purpose GEANT 3.21 package which 
quote as much as 40% discrepancies with ZHS in relevant 
parameters for radio-emission pffl. The main discrep- 
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ancy is on the absolute normalization of the amplitude of 
the electric field, a quantity that is known to be directly 
related to the difference in tracklength between electrons 
and positrons in the shower Jl2[ [l3| . Since the power in 
the signal scales with the square of the electric field am- 
plitude, the discrepancy is unacceptable. The simulation 
of high energy showers in dense media is a cornerstone 
to interpret future measurements and a crucial tool to 
optimize the experiment's performance. Clearly the dis- 
crepancies are numerically very important, they need to 
be understood and the relative merits of alternative cal- 
culations have to be clearly addressed. 

In this paper we compare the ZHS code in the energy 
range between 100 GeV and 10 TeV to simulations per- 
formed with both the GEANT 3.21 package and a more 
recent version of GEANT, namely GEANT4. We find 
GEANT4 predictions to agree with those of ZHS to about 
10%. Moreover, provided that particular care is taken 
in defining some internal variables in the GEANT 3.21 
package, we obtain results consistent with both GEANT 
4 and ZHS. We understand that our result solves the long 
standing discrepancies between the two groups, and gives 
confidence on the calculations, which is essential for the 
future interpretation of data. Our results may be also 
of interest to neutrino experiments which are sensitive to 
tracklengths in water, such as AMANDA (27), BAIKAL 
|§, ANTARES |2{|, SNO M, SuperKamiokande M, 
and even to cosmic ray experiments using water Cerenkov 
tanks such as Auger p|. 

This paper is organized as follows. Firstly we discuss 
the basis of radio interference calculations and remind the 
reader of the physics involved in the coherent Cerenkov 
radio emission from showers. In section III we briefly dis- 
cuss the inputs of the ZHS and GEANT simulators, and 
we quantitatively discuss the differences between the sim- 
ulations. We also explore the main differences between 
the codes, namely the implementation of the LPM ef- 
fect in the GeV- TeV range and its influence on shower 
development. Section [V concludes the paper. 



II. THE PHYSICS OF COHERENT CERENKOV 
RADIO EMISSION AND ITS SIMULATION 

Radio experiments are sensitive to the very low energy 
electrons in the shower (~ 100 keV kinetic energy) which 
contribute most to the excess charge and are responsible 
for the bulk of the coherent radio emission. Hybrid |l^, 
|l(| and analytical methods j2j| have proved very useful 
for calculating the properties of the radio emission, but 
the complexity of the interference phenomena suggests a 
Monte Carlo approach in first instance. 

A realistic simulation of showers for radio applica- 
tions must follow all the particles explicitly down to 
the 100 keV-MeV kinetic energy range. The simulations 
must keep track of the space-time positions of the par- 
ticles in the shower with great accuracy to be able to 
establish the patterns of radio emission at MHz-GHz fre- 



quencies. This is a very challenging problem that was 
first approached by developing a specific package |l^, . 

Once all the shower particles are correctly accounted 
for by the simulations, the contributions to the electric 
field from all the individual tracks with the adequate rel- 
ative phases must be calculated. It is convenient to work 
in the Fraunhofer limit with the Fourier time-transform 
of the electric field amplitude, E(u,x), at a point x in 
the direction of the wavevector k. It has been explic- 
itly derived |l3|, |3j| that the contribution to the electric 
field by a charge e moving at constant velocity v for an 
infinitesimal time interval (ti,ti + <5t) is given by: 
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where eo is the permittivity of the vacuum and fi r the rel- 
ative permeability of the medium, c the speed of light in 
vacuum, lu the angular frequency, and v± the projection 
of the particle's velocity in the direction perpendicular to 
the direction of observation. 

Eq. (|l|) is the basis for all the calculations of the elec- 
tric field amplitude from showers in dense media that 
have been performed so far. Several remarks are worth 
emphasizing. The electric field is proportional to the 
tracklength projected in the direction perpendicular to 
k. Eq. ([j]) has two phase factors. The first phase can 
be written as (1 — n(3 cos 9) ujt\ (n being the refraction 
index and (3 = v/c), and takes care of the relative time 
and position at the onset of the time interval. The sec- 
ond phase is determined by the overall distance between 
the observer and the initial point of the track (R). The 
Cerenkov condition is equivalent to the first phase van- 
ishing and all the points along the track contributing 
coherently to the electric field. 

In order to extend the calculation to many shower 
tracks, the contribution to the electric field from each of 
them as given in Eq. (Q) has to be summed, taking into 
account the different phases which are determined by the 
relative positions of the tracks. To apply this algorithm, 
individual particle tracks in a shower are subdivided in 
sufficiently small intervals. These should be selected so 
that the particle velocity is sufficiently constant for the 
approximation in Eq. (|l]) to be valid. In practice, it has 
been shown that it is possible to consider the average 
velocity and the end points of complete particle tracks 
Jl3| . This reduces the calculation time enormously. In 
Ref. Jl7],[34| it was shown that the approximation under- 
estimates the amplitude at frequencies above 1 GHz but 
it is quite good for lower frequencies. 

Monte Carlo simulations performed in this way render 
a total electric field amplitude which scales accurately 
with the total tracklength in the shower, provided the sec- 
ond phase factor vanishes. This linear behavior holds for 
all angles of observation with respect to the shower axis, 
as long as the frequency is sufficiently low, correspond- 
ing to wavelengths A greater than the dimensions of the 
shower. At the Cerenkov angle the proportionality ex- 
tends to frequencies ~ GHz. Above the GHz range how- 
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ever random phases destroy the scaling behavior. This is 
mainly due to the lateral spread of the shower but there 
are also effects associated with particle trajectories not 
aligned with shower axis and to time delays. 

It is helpful to think of the angular distribution of the 
Cerenkov pulse as the Fourier transform of the longitu- 
dinal development of the excess charge J17J. The main 
diffraction peak appears at the Cerenkov angle. As the 
shower becomes longer, the angular width of the pulse 
becomes smaller. This is particularly important at ener- 
gies above ^PeV in ice, where electromagnetic showers 
suffer the Landau- Pomerancuk-Migdal effect |3j| and be- 
come very elongated. As a consequence the angular dis- 
tribution of the radio pulses becomes very sharply peaked 

As a final remark, Eq. (|l]) corresponds to the Fraun- 
hofer condition which applies when (A vSt) « R 2 . Eq. 
(Q) is clearly valid for any R provided the time interval 
is chosen sufficiently short. It is possible to apply this 
expression to the Fresnel region, one only has to ensure 
that the relative phase factors between tracks are cho- 
sen adequately. This method was used in reference ]ll| 
to calculate the electric field at small distances to the 
shower. 



III. SIMULATION OF ELECTROMAGNETIC 
SHOWERS: ZHS VS GEANT 

A reliable Monte Carlo (MC) code that simulates 
shower particles in detail is an essential intermediate step 
to estimate the expected radio-emission from a high en- 
ergy _showc_r._Sevcral options can be used for this purpose 
, B9fl and here we concentrate on two of 
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them namely, GEANT in the two versions: GEANT 3.21 
H and GEANT4 |§ and the ZHS code Jl|. We briefly 
discuss the most relevant processes considered by both. 
Further details can be obtained from the references. 

The ZHS is a fast code optimized for the simulation 
of electromagnetic showers in homogeneous ice. It fol- 
lows electrons, positrons, and photons until an exter- 
nally fixed energy threshold is reached. The program 
accounts for brcmsstrahlung and pair production, tak- 
ing screening corrections into consideration as well as the 
Landau-Pomerancuck-Migdal (LPM) effect ]35j . It also 
accounts for multiple elastic scattering, Moller, Bhabha, 
and Compton scattering, as well as electron-positron an- 
nihilation. The time delays of the particles with respect 
to an ideal particle moving along the shower axis at the 
speed of light are carefully implemented, accounting for 
subluminal velocities, geometrical effects, and multiple 
elastic scattering. The longitudinal profile of electromag- 
netic showers obtained with ZHS was shown to agree with 
standard parameterizations of the depth development in 
p3| . In H [l2| ZHS was also shown to agree with standard 
parameterizations of the lateral distribution of electrons 
and positrons in electromagnetic showers. 

GEANT is a well known, well tested and widely used 



simulation and detection package which is suitable for 
a large range of applications. Version 4 is a recent up- 
date with several improvements. Many computational 
parameters such as energy thresholds, or the minimum 
step length in particle propagation, and even cross sec- 
tions and corrections may be included or not by the user. 
However, care has to be taken in setting input parame- 
ters which are consistent and in selecting the appropiate 
physical processes to the problem in consideration p8[ . 



A. ZHS vs GEANT3: Discrepancies 

The most important quantities in the calculation of the 
coherent radiopulses are the total, total projected, and 

The total tracklength is defined 



excess tracklengths 13 



as the sum of the lengths of all charged particle trajec- 
tories in the shower. The total projected tracklength is 
the sum of the projections of the particle paths along 
the shower axis. The excess charge track length is sim- 
ply the difference between positive and negative charged 
projected tracklengths. They are calculated in the ZHS 
code by default. These magnitudes have been found to 
be linearly dependent on the primary energy to a great 
degree of accuracy, and are directly proportional to the 
peak value of the electric field at the Cerenkov angle. 
This is just a consequence of energy conservation in the 
shower. In fact it has been shown that shower fluctu- 
ations do not affect the total tracklength and that by 
looking at the peak value of the coherent electric field, 
an accurate measurement of the shower energy can be 
obtained H |l§ |l§ @ . 

Simulations carried out with GEANT 3.21 in @ pre- 
dict a total, total projected and excess projected track- 
lengths which are between 20% and 40% smaller than the 
corresponding values obtained with ZHS. The number of 
particles at shower maximum in Q is about 35% smaller 
than for the ZHS. In the following we address the origin 
of the discrepancies and establish that the simulations in 
p6[ underestimate the value of these crucial parameters. 

We have simulated showers of different energies us- 
ing ZHS, GEANT 3.21 and GEANT 4 and studied their 
depth development, lateral distribution and the three 
tracklengths which are relevant for radio calculations. In 
order to understand the results of |2(| we have made 
two different sets of simulations changing an internal pa- 
rameter of GEANT 3.21. The parameter is called IA- 
BAN and it is discussed in more detail in the appendix. 
Here it should be sufficient to say that changing the value 
of IABAN is equivalent to fixing different values of the 
external energy threshold below which electrons are no 
longer followed in the simulations. This, we believe, is the 
ultimate cause of the discrepancies in the tracklength be- 
tween ZHS and the results presented in |2(| . We perform 
two sets of simulations with GEANT 3.21: one corre- 
sponding to the default values, set (a), and the other 
using a modified setting which correctly accounts for all 
the electrons in the shower which are above the kinetic 
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energy threshold, set (b), (see the appendix for details). 

Table [| summarizes the numerical values of the total 
tracklength, total projected tracklength, and excess pro- 
jected tracklength in ice as obtained in GEANT 3.21, 
GEANT 4, and ZHS. We averaged the tracklength over 
100 electron-induced showers of energy 100 GeV. The 
electron kinetic energy threshold is -£Tth=100 keV (or 
Bth = 0.611 MeV total energy) in all simulations. Us- 
ing the default version of GEANT, GEANT 3.21 (a), 
we reproduce the results of [^6) to the statistical ac- 
curacy of the simulations. With the setting GEANT 
3.21 (b), which accounts for an external threshold Eth = 
0.611 MeV consistent with what is explicitly done in the 
ZHS MC, the 40% discrepancies between [g6| and ZHS 
on the tracklength as well as on the number of parti- 
cles at shower maximum, are reduced to 10%. Further- 
more the results of ZHS, GEANT3 (b) and GEANT 4 
are within 10% of each other. In Ref. |2(| it is argued 
that an upper bound to the tracklength per unit en- 
ergy is simply obtained by calculating the inverse of the 
minimum energy loss of an electron in ice [dE/dX] m l n . 
We remark here that this is indeed an upper bound to 
the average length traveled by a single electron, how- 
ever it does not take into consideration the contributions 
to shower tracklength through secondary electrons pro- 
duced by bremsstrahlung photons that contribute to the 
energy loss of the primary electron. 
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FIG. 1: Average number of electrons plus positrons in 100 
GeV electron initiated showers, averaged over 100 showers, as 
a function of depth in ice. The profiles were obtained with dif- 
ferent Monte Carlo programs; ZHS (squares), GEANT 4 (un- 
filled triangles), GEANT 3 with IABAN=2 (filled triangles), 
GEANT 3 with IABAN=1 (filled inverse triangles). Depth is 
measured in radiation lengths ( 1 r.l.= 36.08 g/cm 2 ). 

Fig. [l] (Fig. ||) shows the average longitudinal devel- 
opment of the total (excess) charge for 100 electron- 
initiated showers of energy 100 GeV in ice. The same 
MC codes as in table | were used to produce the plots. 
The number of particles at shower maximum as obtained 
with GEANT 3.21 (a) (see table caption and appendix) 
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FIG. 2: Same as Fig. jj] for the average number of electrons 
minus positrons. 



agrees with the result quoted in ]26| , and differs by 30% 
from the value predicted by ZHS. However the expec- 
tations from GEANT 3.21 (b), GEANT4, and ZHS are 
again within ~ 10% of each other. In Figs. ||, [| and || we 
show the lateral distribution of electrons and positrons at 
shower maximum corresponding to ~ 7 radiation lengths 
(1 r.l. = 36.08 g/cm 2 ), as well as at 4 r.L, and at 10 
r.l. for the four nominal MC codes. Again, GEANT 4, 
ZHS, and GEANT 3.21 (b) give similar results whereas 
GEANT 3.21 (a) underestimates the total number of par- 
ticles and also slightly distorts the lateral distribution. 
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FIG. 3: Lateral distribution of electrons and positrons in 100 
GeV showers in ice at a depth of 7 radiation lengths as ob- 
tained with different MC codes. ZHS (squares), GEANT 4 
(unfilled triangles), GEANT 3 with IABAN=2 (filled trian- 
gles), GEANT 3 with IABAN=1 (filled inverse triangles). All 
densities are normalized to the same number of particles. 1 
Moliere unit is 10.4 g/cm 2 in ice. 



The discrepancies between the radio emission Simula- 
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TABLE I: Total, total projected and excess projected tracklength in ice as obtained in GEANT 3.21, GEANT 4 and ZHS. The 
results show the average over 100 electron-induced showers of primary energy 100 GeV. The electron kinetic energy threshold 
is 100 keV (0.611 MeV total energy) in all simulations. GEANT 3.21 simulations were performed with the following values of 
IABAN (see appendix): (a) IABAN=1 and (b) IABAN=2. The numbers in parenthesis indicate the relative differences (in %) 
taking the GEANT4 results as reference. 



MC Code 


GEANT3 (a) 


GEANT3 (b) 


GEANT4 


ZHS 


Total track [m] 


409.2 (-30) 


577.9 (-2) 


587.9 (0) 


642.3 (9) 


Total projected [m] 


372.9 (-18) 


450.0 (-1) 


453.2 (0) 


516.7 (14) 


Excess projected [m] 


92.4 (-25) 


123.5 (-1) 


122.7 (0) 


132.4 (8) 


Excess/Total 


0.225 (+8) 


0.214 (+2) 


0.208 (0) 


0.206 (-1) 
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FIG. 4: Same as figure W at a depth of 4 radiation lengths. 
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overall behavior obtained by ZHS up to 5 GHz is very 
similar to that shown in ref . p6[ , provided the normal- 
ization is corrected according to the tracklength deficit. 
It is however difficult to be confident on the results at 
higher frequencies, because very small relative time dif- 
ferences of only 0.02 ns and distances between particles 
of the order of 6 mm, are expected to have an important 
effect on the interference patterns. Theoretical calcula- 
tions based on simple current density parameterizations 
which neglect time delays and front curvature Jl8|, [2^] are 
not expected to be useful in this frequency range for the 
same reason. 

The good agreement between GEANT 3.21 (b), 
GEANT 4 and ZHS is a remarkable achievement given 
the fact that ZHS and GEANT have been developed in- 
dependently of each other and are of very different na- 
ture and purpose. There are two important differences 
between the codes but a detailed analysis is not the pur- 
pose of this paper. In the next subsection we show that 
none of them can account for a 50% discrepancy in the 
tracklength, although they might explain the remaining 
(~ 10%) discrepancies in the observables. 

Finally we would like to stress that the results pre- 
sented above are all reproducible with the available (upon 
request) codes [[S0|. We believe we have solved this long 
standing and crucial discrepancy between the different 
simulations. The implications for radio projects are ob- 
vious but it may be stressed that other simulations re- 
lying on tracklength in water detectors that have been 
made for many other experiments such as those com- 
puting Cerenkov light in water/ice, should be checked 
against the results presented in this paper. 



B. ZHS vs GEANT: Differences 



FIG. 5: Same as figure W at a depth of 10 radiation lengths. 



tions presented in our previous work jl2|, [l^, [l5| [l6|, [Tt], 
|l8f and those reported in |^6) can be mostly interpreted 
in terms of the normalization of the tracklength. The 
remaining differences can be attributed to the effect of 
track subdivision. If small subintervals are chosen, the 



It should be remarked that both codes take the same 
interactions into account, except for two important dif- 
ferences, namely the LPM corrections which are not ac- 
counted for in GEANT 3.21 (39), and the photoelectric 
effect which is ignored in the ZHS code. The photoelec- 
tric effect is not expected to be important for low Z nu- 
clei such as hydrogen and oxygen for energies above the 
electron Cerenkov energy threshold (~ 100 keV kinetic 
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energy in ice). We will discuss the LPM effect and show 
that it hardly represents any effect for the tracklength 
calculations. 

The LPM effect stems from the fact that the inter- 
action distance in bremsstrahlung or pair production is 
inversely related to the momentum transferred to the 
nucleus (q). When the initial and final electron mo- 
menta in bremsstrahlung (or the electron and positron 
momenta in pair production) become ultrarelativistic, 
the two electron momenta are almost collinear and the 
longitudinal momentum transfer qu can be very small. 
Conversely, the distance along which the interaction oc- 
curs might become large, comparable to the interatomic 
spacing which depends on the density of the medium. 
Under these circumstances the electron encounters addi- 
tional atoms along the interaction distance which cause 
multiple Coulomb scattering. This introduces destruc- 
tive interference in the interaction matrix element and 
as a result the total bremsstrahlung and pair produc- 
tion cross sections are suppressed Q. The suppression 
becomes important when the average multiple Coulomb 
scattering angle (9 S ) is comparable to the scattering an- 
gle in the bremsstrahlung or pair production interactions 

In pair production the two final electron momenta are 
collinear only when the photon energy is above an en- 
ergy scale (£"lpm) which decreases with the density of 
the medium (Elpm ~ 2 PeV in ice). The same occurs in 
bremsstrahlung when the energy of the electron is above 
Elpm, but also when the fraction of the electron's energy 
carried away by the photon is small. 

Following Ref. [Q a parameter (s) is introduced 
through the condition 9 S ~ 6*;. When s < 1 the LPM 
has an important effect on bremsstrahlung and pair pro- 
duction interactions, s is given by pO): 



s v \ I — -ir— bremsstrahlung, (2) 



E 1-U 



s oc 



' EhPM 1 



E v(l-v) 



pair production, (3) 



where u is the fraction of the energy of the electron car- 
ried by the radiated photon in bremsstrahlung, and v is 
the fraction of the photon's energy carried by the elec- 
tron or the positron in pair production. Clearly s < 1 
in both equations when E 3> Slpm but also when u is 
smaller than mlpm — E/Elpm in bremsstrahlung. 

As a consequence, it can be easily shown that the dif- 
ferential bremsstrahlung cross section scales as dN/dk oc 
l/Vk instead of exhibiting the Bethe-Heitler behavior 
dN/dk oc where k is the energy of the photon. The 
integral of the differential cross section is naturally reg- 
ularized, preventing the infrared catastrophe at very low 
photon energies. Moreover, the average fraction of the 
electron's energy carried away by bremsstrahlung pho- 
tons (it) increases with energy due to the suppression of 
photons with energy below ulpm-E. This effect has been 



experimentally measured at energies ~ 10 — 25 GeV in a 
variety of materials |4^, Q (see 0] for a comprehensive 
review) . 
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FIG. 6: Total bremsstrahlung cross section as a function of 
electron energy. The solid line is the total cross section when 
the LPM effect is accounted for. The dashed line is the Bethe- 
Heitler cross section i.e. without the LPM effect. The dashed- 
dotted line is the total cross section as obtained integrating 
the differential cross section in reference ^] . The cutoff pho- 
ton energy is 100 keV for the three curves. 

Fig. [6] shows the bremsstrahlung cross section in ice 
when the LPM effect is considered and the well known 
Bethe-Heitler cross section (without LPM effects). De- 
spite the fact that -Elpm is in the PeV energy range the 
suppression of the bremsstrahlung cross section starts 
already at around 10 GeV due to the effect explained 
above. The LPM bremsstrahlung cross section drops as 
y~E in contrast to the logarithmic increase with energy 
of the Bethe-Heitler cross section. Also shown for com- 
parison is the LPM cross section obtained in p0[ ]. The 
small differences between both LPM cross sections can 
be attributed to the Koch and Motz empirical low en- 
ergy corrections Ej considered in ZHS but ignored in 



Shower development is nevertheless fairly insensitive 
to the large differences between the LPM and the Bethe- 
Heitler bremsstrahlung cross section in the energy range 
below 100 TeV. We have calculated the electron energy 
loss per radiation length in the 100 GeV - 100 TeV energy 
and we have seen that it is unaffected by the inclusion or 
not of the LPM effect. The reason for this is that the de- 
crease in the bremsstrahlung cross section in this energy 
range, is compensated by the expected increase in (u), so 
that the product (u)a (i.e. the energy loss) is the same 
as in the Bethe-Heitler calculation. In contrast, at en- 
ergies above -Elpm this cancellation does no longer exist 
and the energy loss including the LPM effect is lower than 
the one predicted by the Bethe-Heitler cross section. This 
added to the fact that the pair production cross section is 
also suppressed above -Elpm, leads to the characteristic 
elongated longitudinal profiles of LPM electromagnetic 
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showers 15, 46 . 



Finally it is worth remarking that the total tracklcngth 
is completely insensitive to the LPM effect at all energies. 
This is just a manifestation of energy conservation in the 
shower as explained before. 

There are some other minor differences in the treat- 
ments and approximations used for the cross sections for 
relevant electromagnetic processes which have a slightly 
different implementation in GEANT and ZHS. The treat- 
ment of Moliere multiple elastic scattering is also differ- 
ent. However in Ref. |26| it is shown that a small varia- 
tion in these cross sections by as much as 25% cannot be 
responsible for large discrepancies in shower parameters. 
Surely these minor differences together with the LPM ef- 
fect are likely to be responsible for the remaining 10% 
differences between ZHS, GEANT 3.21 (b) and GEANT 
4. 



Small discrepancies that remain between ZHS and 
GEANT are possibly due to minor differences which can 
be attributed to physical uncertainties in the processes 
involved. As a result all previously published calcula- 
tions based on the ZHS Monte Carlo are expected to be 
correct to the same degree of accuracy, modulo other un- 
certainties which were discussed in the original papers. 
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IV. CONCLUSIONS 



We have simulated electromagnetic showers in ice using 
two different Monte Carlo codes, GEANT and ZHS. The 
discrepancies discussed in reference 1 26 between the ZHS 
code and GEANT diminish to the 10 % level when the 
GEANT 4 version is used, and also when the GEANT 
3.21 code is carefully set to allow for a constant energy 
threshold. 

We have shown that the influence of the LPM effect 
on shower development is negligible for energies below 
E-lpm, despite the fact that the bremsstrahlung cross sec- 
tion is significantly affected. As it has been now well es- 
tablished, the shower longitudinal profile is dramatically 
affected by the LPM at energies above £Lpm- 

It is worth stressing that ZHS is a Monte Carlo specif- 
ically designed to calculate coherent radiopulses in elec- 
tromagnetic showers. On the other hand, GEANT is 
a general purpose Monte Carlo designed to deal with a 
large variety of problems. The fact that both unrelated 
and independent programs give results to the 10% level 
of precision gives us confidence on both packages. The 
ZHS program is also suitable for simulations of experi- 
ments that use ice or water and the Cerenkov imaging 
technique, as long as muons are not important. It can 
also be extended to other media. In addition, ZHS is 
considerably faster than the GEANT code and can reach 
much higher energies, making it adequate for high energy 
neutrino experiments. 



APPENDIX A: THE GEANT SETTING IABAN 

The default setting of GEANT includes a variable 
named IABAN |3l] whose value is set to 1. The effect 
of this setting is to stop tracking a fraction of the elec- 
trons and positrons of energy below 10 MeV, even when 
the external electron energy threshold is set to a lower 
value. To be more precise, when the distance to the next 
bremsstrahlung interaction is larger than the range of 
the electron, the electron is stopped and its tracklength 
is not computed. Typically in a 100 GeV shower ~ 50% 
of the electrons of energy below 10 MeV are stopped and 
eliminated from the simulation. The stopped electrons 
account for 40% of the total tracklength. As a conse- 
quence the total, total projected and excess projected 
tracklengths and hence the total electric field emitted by 
a shower are underestimated by ~ 40%. 

If the variable IABAN is set to 2, all the electrons and 
positrons in a GEANT simulation are tracked down to 
the external energy threshold which is chosen by the user. 
This setting gives the correct tracklength. 

In GEANT4 the switch IABAN does no longer exist, 
and particles are always correctly tracked down to the 
external threshold. As a consequence all the particles 
that build up the radio-emission are correctly accounted 
for. 
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